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In this article, we study control of electron-nuclear spin dynamics at magnetic field strengths 
where the Larmor frequency of the nucleus is comparable to the hyperfine coupling strength. The 
quantization axis for the nuclear spin differs from the static Bq field direction and depends on the 
state of the electron spin. The quantization axis can be switched by flipping the state of electron 
spin, allowing for universal control on nuclear spin states. We show that by performing a sequence 
. . . of flips (each followed by a suitable delay), we can perform any desired rotation on the nuclear 

' spins, which can also be conditioned on the state of the electron spin. These operations, combined 

, with electron spin rotations can be used to synthesize any unitary transformation on the coupled 

■ electron-nuclear spin system. We discuss how these methods can be used for design of experiments 

for transfer of polarization from the electron to the nuclear spins. 

PACS numbers: 03.67.Lx 
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I. INTRODUCTION 

In this paper, we study the problem of control of coupled electron-nuclear spin system consisting of an electron, 
coupled with one or more nuclear spins. Manipulation of such electron-nuclear spin systems is fundamental in the 
field of ESR/EPR 1] with application to the study of structure and dynamics of paramagnetic species. Coupled 
electron- nuclear spin systems have recently been studied from the perspective of quantum information processing 

Study of control of these systems is interesting from perspective of 



, quantum control as dynamics of coupled electron-nuclear system has some salient differences compared to coupled 



nuclear spin systems 



17|. 



. The Rabi frequency of the electron (at typically available microwave power in pulsed EPR experiments) is of orders 
of magnitude larger than the Rabi frequency of the nucleus (at typical rf-power) . Duration of pulses on the electron 
- ■— , , and nucleus are in ns and /is regimes respectively in EPR experiments. Local rotations on electrons, therefore take 
much less time than the nucleus. The energy eigenstates of the nucleus (at field strengths where hyperfine interactions 
and Larmor frequency of the spins is comparable) depend on the state of the electron and are not aligned with the 
static i?o field. This results in the well studied phenomenon of Electron Spin Echo Envelop Modulation (ESEEM) 
This opens the possibility of controlling the nucleus by simply manipulating the electron spin between its two 
eigenstates (by a series of vr pulses and delays). The quantization axis and the precession frequency (around this axis) 
of the nucleus is switched every time the state of electron spin is flipped, making the nucleus nutate around the new 
quantization axis. We will shown that by performing a sequence of flips on the electron (each followed by a suitable 
delay), we can perform any desired rotation on the nuclear spin states. This rotation can also be conditioned on the 
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state of the electron spin. This mode of control on nuclear spins obtained by switching between two rotation axis 
is an excellent example of so called switched control systems in control theory [l^ . Controlling nuclei by switching 
quantization axis is preferable than rotating nuclei with rf-fields as hyperfine couplings are orders of magnitude larger 
that Rabi frequency obtainable with typical radio frequency power. Therefore much shorter experiments can be 
designed, which can significantly reduce relaxation losses, mainly arising from short transverse relaxation times of 
electron. Furthermore, since the only manipulation required involves flipping the state of the electron spin, we show 
that it is easy to design these experiments so that they are robust to the Larmor dispersion of the electrons. It is 
expected that this switched mode control of nuclear spins will find applications in quantum information processing 
and pulsed EPR, including pulsed dynamic nuclear polarization (DNP) experiments . 

II. THEORY 



We consider as our model, a spin system consisting of one electron spin S — ^ and one nuclear spin I — (see, e.g.. 
Sect. 3.5 of [1]) in a static magnetic field Bq along the z direction. We will later also discuss the case of one electron 
coupled to many nuclear spins. The Hamiltonian Hq of the electron- nuclear spin system in the laboratory frame may 
be written as 



Hq = UJsSz + iJJlh + S • A • I, 

where los represents the Larmor precession frequency of the electron spin, uji represents Larmor precession frequency 
of the nuclear spin, and A is the electron-nuclear hyperfine coupling tensor. Given the Pauli matrices Ux '■— ( i o)i 
ay i (i V) ' ^-'^^ '■~ (o -i) ^i^d the identity matrix cto := (J i)j the operators Sj and Ik {j,k G {x,y,z}) are 
defined by Sj — [dj ®(Jq)/2 and Ik = (cto ^ o'fe)/2 (see |l7ji). The static field i?o, sets the quantization axis of electron 
spin and the coupling Hamiltonian = S ■ A ■ I, averages to 

He — ASzIz + BxSzIx + BySzIy, 

The transverse plane axis for the nuclear spin subspace can be so chosen that the last two terms are combined to 
form 



H, = ASJz+BSzIx, 



where B = ^ B'^ + B'^. The full Hamiltonian of the system then takes the form 



Hq = UJsSz + i^lh + ASzh + BSzIx- 
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Here ojs = —JcBq and uji = — 7„i?o- Here 7e is the gyromagnetic ratio of the electron (negative) and 7„ is the 
gyromagnetic ratio of the nucleus (we take positive, as for a proton). 

Let a and /3 denote the state of the electron spin oriented along and opposite to the direction of the static magnetic 
field Bq respectively. For an electron, the a configuration has higher energy than the (3 configuration. 

When the electron is in the a state, the nucleus sees a net field 



A B ^ 

When the electron is in the j3 state, the nucleus sees a net field 

/ A , B 

Bi3 = {Bo + —)z + —x 
27« 27„ 

The nuclear precession frequency in the two states is given by 
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(1) 



where ajj is negative. Let da and df3 denote unit vectors along the field directions Ba and Bp. The states \aa) and 
\a(3) represent the states of nuclei oriented along or against da respectively, when the electron is in the a state. These 
states constitute the a manifold. Similarly, the states \(3a) and \f3(3) represent the states of nuclei oriented along or 
against the field dp respectively, when the electron is in the /3 state. These states constitute the /? manifold. Fig. [T] 



depicts the field da and dp and the energy eigenstates states when > 
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FIG. 1: Fig. A depicts the vectors da and rf^ when \uji\ > ^. Fig. B depicts the corresponding energy level diagram. 



In the rotating frame (rotating with electron at frequency Wg), the Hamiltonian of the electron- nuclear spin system 
takes the form 



H, = nsS, + ujih + AS J, + BS.h 
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where Jig is the resonance offset for the electron. We assume for now that precession frequency of the electron is 
well defined and ilg — (the case when fig will be discussed subsequently). We now show how by switching 
between the a and (3 states of the electron, by tt pulses, we can synthesize any desired rotation on the spin states 
of the nuclei. Furthermore, since uja ^ w^, we will show these rotations can be made conditioned on the electronic 
state. To fix ideas, consider the state transformation 



\(3a) ^ m- 



When the electron is flipped from f3 to a state, the effective field felt by the nucleus switches from direction dp to 
da and the nuclei begin to precess around this new axis. Following this precession for time ti, we flip the electron 
back, causing the precession axis to return to dp and then let the precession happen for time T2 and so on. Since da 
and dp are two independent axis of rotations, switching between them followed by precession can generate any three 
dimension rotation of the nuclear orientation. Fig. shows the trajectories of nuclear spin on the Bloch sphere as 
it precesses around da and dp directions. These trajectory plots can be used to construct the sequence of switchings 
between the axis (and delays between switchings) da and dp to produce the desired rotation of nuclear spin state. 
Fig. [2j3, shows one such graphical construction for finding the switching times ti, T2 and T3, that rotates the unit 
vector da to vector n, by alternate rotations around dp and da- 

dn 





B 



FIG. 2: Fig. [2}V shows the various trajectories of nuclear spin on the Bloch sphere as it precesses around da and dp directions. 
This trajectory plot aids in visual construction of rotation between points on the sphere by alternation between rotations 
around da and d/s. Backtracking these trajectories can help to construct desired rotation of nuclear spin state by switching 
between field directions da and d/s. Fig. [2)3 depicts one such construction (projected on the plane of axis da and d/s ) that 
rotates unit vector da to the unit direction n. The initial rotation is around axis dfs for time n, followed by rotation around 
da for time T2 and finally around dp for time ra. The times ri, T2 and T3 can be explicitly computed from this figure. 
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Fig. [3] shows the transformation \f3a) with a three period pulse sequence as shown in the Fig. [3l The 

delays ti, T2 and T3 are chosen so that the nuclear spin magnetization execute the trajectory shown in Fig. The 
last TT pulse in the figure returns the electron to state p. The times ri, T2 and can be readily computed from this 
figure. 




B 
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FIG. 3: The Fig. A depicts the trajectory executed by nuclear spin magnetization, initially oriented along the direction d'', 
when the pulse sequence in Fig. B is implemented. The transformation is the inversion of the nuclear spin resulting in the 
transformation \/3a) —> |/3/3). Fig. C, shows the desired rotation (shown with solid arrow) \/3a) —> is implemented by 

switching the electron between a and /3 manifold. 

Similar arguments can be used to transform any energy eigenstate into another by simply flipping the electron 
spins followed by appropriate delays. This then naturally leads us to the question of number of times the state of 
the electron needs to be switched to produce any desired rotation on the nuclear spins, in a given eigenstate in either 
a or (3 manifold. The problem of flnding the minimum number of switchings for producing arbitrary rotations by 
switching between two axis of rotation has been studied before. We recapitulate the main results and motivate further 
extcntions. 

Remark 1 Uniform Generation: Let 7 denote the angle between two rotation axis da and dp as shown in Fig. ([3]). 
Then it is a known result 20] that any three dimensional rotation in group S0{3) (similar result holds for group 
SU{2)) can be constructed by not more that k + 2 rotations (performed around da or dp), where 

TT TT 
k+1 - ^ k 

Therefore, starting in cither a or (3 manifold , any rotation on nuclear spin can be obtained by a pulse sequence as 
shown in Fig.(I3j3) with no more that k + 2 evolution periods Ti. Furthermore, there exist algorithms to find the 



sequence with minimum number of switches [21| to generate any arbitrary rotation using only rotations around da 
and d^. Here, we present a strategy for obtaining such a sequence (which need not be minimal in terms of number of 
switches) which is easy to construct geometrically. 

Let Ti — Tfe (k even) be the times for alternate rotations (starting with c?^) around and da, with angular frequency 
ujf^ and LOa, so chosen that the resulting rotation takes the initial vector da to unit direction n as in figure (I2j3). Let 
C/(a,T) = exp(— iwoT da ■ a) and U{I3,t) — exp(— iw/jr dp ■ a), be rotations in SU{2) around axis da and dp. A 
rotation by angle 6 around an unit vector h — rixi + Uyj + n^fc, takes the form 

9 

exp{-'i9n -cr) = cos— 1 — i 2 sin — n • cr, (2) 
where n ■ a = UxUx + nyUy + nz<Jz- Then, a sequence of rotations 

C/(a, Tk)... U{p, n) Uia, t) U{(3, -n) . . . C/(a, -Tfc), (3) 

V ' " v ' 

will generate a rotation around the unit vector fi with an angle 9 — TLJa- Let 

27r 27r 

ta = — ; tfi = — ; (4) 
The rotations U{(3, — r) and U{a, —t) are simply obtained by observing that 

[/(/?, i^-r) = -[/(/?, -r); U{a,ta - t) = -U{a,-T) (5) 

Since k is even, the negative signs in ([3]) cancel out (if k is odd, it just produces a global phase). We now show how 
all these results need to be further generalized if one is interested in producing controlled operations on the nuclei. 

A. Controlled Operations on Nuclear Spins 

A a manifold selective controlled operation on the nuclear spins transforms \aa) — > \aa') and jo;/?) — > |a/3'). While 
the (3 electron manifold is left unperturbed, i.e., \(3a) \(3a) and |/3/3) — ^ \Pf3). We show that such controlled 
operations can again be performed by toggling the state of the electron spin with tt pulses, followed by suitable delays. 
This is best seen by observing that by flipping the state of the electron spins, we can switch between the effective 
Hamiltonians (generators) 



H^=ujih+SMlz+BIx) 



(6) 



and 



(7) 



For example, the pulse sequence ti — tt — T2 — tt (where tt pulse is on the electron) evolves Hi for time ti followed 
by H2 for time T2. A straightforward calculation shows that when LUa ^ top (in Eq. [T]), the repeated commutators of 
—iHi and —iH2 generate a six dimensional algebra, spanned by generators t — — ly, Iz, Szlz, Szlx, Szly}. Then 



a standard result on controllability 



states that any unitary transformation generated by these generators can be 



produced by simply switching between —iHi and —iH2, i.e., any rotation of the kind 

eyi^{~i{aSzIp + bIq)), (8) 

where Ip and Iq are nuclear spin operators can be synthesized. Specifically, we can perform a rotation exp(— zS'"/-y) 
that only rotates nuclei in the electron's a manifold {S" — (^ + Sz)^ S'^ = (i — Sz) )■ 
We rewrite, 



Hi = ujiiSf + Sl)Iz + ''-^^^^{AIz + BI,), 

which when rewritten looks like 



which written in a 6 x 6 matrix form takes the block diagonal form 



„ Do. y rr - 



where denotes rotation around axis with angular frequency cJq, and Dp denotes rotation around axis dp with 
angular frequency LOp. The algebra t is simply all block diagonal 6x6 skew Hermitian matrices, i.e., matrices of the 
form \ ? ) , i.e. Ai ~ —A\ and A2 ~ -~A\. Note, in the special case when A = in ([T]), then LOa = LUp in 



In this special case the algebra generated by Hi and H2 is only 3 dimesional algebra and we cannot produce arbitrary 
controlled operations on the nuclei. This case however is of only theoretical interest and we will not consider it any 
further. 

A TT pulse on the electron will switch Hi to H2 and vice- versa. A sequence of delays and flips ri — tt — T2 — • • • — Tfe — tt, 
where k is say even, will execute the following rotations on the nuclear spins in a and /3 manifold respectively 



Uc. = U{(3, Tfc) . . . U{a, Ti); Up = U{a, r^) . . . U{f3, n). 

Here U{a,T) = exp(— iZJ^r) and U{/3,t) — exp{~iDpT). The goal is to find the sequence of times tj,, such that 
Up — I and Ua — Uf is the desired rotation. This constitutes a controlled rotation on nuclear spins, conditioned on 
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the electron state. Similarly, finding a sequence of times r/c, such that Ufs = Ua = Uf constitutes a uniform rotation 
on the nuclei. In general, we want to find a sequence that performs manifold selective operations, i.e., Ua and Up are 
as desired. It is also a known result [20] that there exists a finite number of switchings N, such that arbitrary Ua and 
Ufj can be constructed by no more than N number of switchings. Finding methods to generate these in minimum 
number of switchings is an interesting problem. We present a constructive search strategy for achieving this (our 
methods need not be minimal in terms of number of switching) . 

Let us perform a series of rotations on the (3 manifold (with k even) given by 

Up = f/(/3, -n) . . . Uia, -Tfc) Uia, r^) . . . i7(/3, n), 

V V ' 

where t/g is implemented as 

Up = U{P, tp-n)... U{a, - Tfc) U{a, n-) . . . i7(/3, n). 

By definition. Up = \ . We now have the independence in selecting ti — t^, in such a way that the corresponding 
rotation on the a manifold, given by 

U{a,tp-Ti)...U{P,t^-Tk) Ui(3,Tk)...U{a,Ti) =Uf. (9) 

We rewrite this relation as 

1 = UiP, Tk-U)... Uia, Ti - tp) Uf U{a, -n) . . . C/(/3, ~n), (10) 
^ ^ ^ ^ ^ ^ 

Equation (fTO|) . can then be written as, 

1 - BU{f3, Tfe) . . . AUia, Ti) Uf U{a, -n) C/(/3, -n), (11) 



where A = U(a, —tp) and B ~ U{(3, —ta). 

Equation (jlip can be used to determine ti — Tfe. In equation (llip . the times ti — r^, are chosen to increase the trace 
of Up to 2 (the only matrix in SU{2) with trace 2 is 1 ). These can be obtained iteratively in a closed form by the 
following optimization procedure. Suppose T2 — Tk is fixed. Now the Eq. PH) . can be written as 

rr(C/|) = Tr{C Uia, n) Uf U{a, -n)), 

where the time ti is chosen to maximize Tr{Up). This optimization then takes the form 



tr{CU{a,Ti)UFU{a,-Ti)) = 2{cos - cos y -sin-siny7ic •'^(ti)}, (12) 



where Uf = exp(— i0n ■ cr) and C = exp(— i^c'^c • cr). 
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In the above expression ti generates rotation around axis da and only rotates n±, the part of n, that is perpendicular 
to da- The above expression is maximized when is aligned or anti- aligned (depending on the sign of sin | sin ^) 
with {nc)_L, the part of (nc)± that is perpendicular to d^. We can write an explicit expression for ti. Let a be a 
unit vector in the direction of sin |(n — (n • da)da), and 6 be a unit vector in the direction of sin ^{nc — {ric ■ da)da), 
then sin(Ti) = (b x a) ■ da- Having determined the optimal value ti, we can now proceed to maximize T2 which 
is now an optimization of a function like Tr(Ci U{I3,T2) Up U{j3, —T2)), and so on. Having found a set of values 
Ti — Tfe, this way, we iterate again. Each iteration increases the value of tr{Up), and therefore this value stabilizes to 
a limit. If this value is not 2, we can increase k and repeat the iteration. The resulting set of evolution times ti — 
synthesize a selective rotation on the a manifold as in Eq. ([9]). Similarly we can synthesize a selective rotation on 
the (3 manifold and therefore any rotation of the kind given in Eq. ([5]). Alternatively, we can synthesize a desired 
unitary transformation U(} on the beta manifold as described in Eq. ([3]), which results in some rotation on the 
alpha manifold (that is easily computed). Now performing an a manifold selective rotation of UaiU^)^^ , will result 
in the the desired rotations on the two manifolds. Using the above described method, we explicitly compute the pulse 
sequence for implementing a controlled rotation exp(— i7rS'"/z), which inverts the nuclear spins in the a manifold and 
leaves the (3 manifold unaffected. We assume ^ — .5 and ^ — 1. This gives — 'in and tp = The pulse 
sequence 

Tl - TT - T2 . . . T5 - TT - ta - TT - {t/S - T5) - TT - {ta - T4) ~ . . . ~ TT ~ {tp - Tl), 
^ V ' 

P 

where the switching times, in the units of 1/A are ti = 10.32, T2 = 4.64, = 9.75, T4 = .2760, T5 = 10.97 and 
Tg = 5.47, implements the desired rotation. Also note that using the above sequence as a building block, the pulse 
sequence P — tt — P — tt implements the rotation exp(— ztt/^). Similarly we can find a switching sequence Q, that 
implements exp(— i7r/^). Then the pulse sequence P — Q — P — Q will implement the rotation eyi'p{— 11:21 z) ■ 

Remark 2 (Synthesizing Arbitrary Unitary Transformations:) Any unitary transformation on the electron- 
nuclear spin system is given by |23| 

UiUs eM-i[aS.J^ + bSyly + cSJ,])ViVs 

where Ui, Vj are single spin operations on the nucleus and Us, Vs are single spin operations on the electron. The 
operators, IxSx, lySy, IzSz all commute. 

As described in Eq. ([8]), we have shown that all rotations of the kind eyLY>{—iOiSzIa) and exp(— i02^h) can be 
synthesized by simply toggling electron spins (here la and are arbitrary nuclear spin operators). These rotations, 
along with single spin rotations, exp(— z^aS-y) on the electron can be used to perform arbitrary unitary rotation on 
the electron-nuclear spin system. The bilinear rotations exp{~i9SxIx) can be synthesized as 

exp{-i6SxIx) = exp{~i^Sy)exp{-i9SzIx)exp{i^Sy). 
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Remark 3 (Refocusing Resonance Offsets for Electron Spins:) In the above described methods, we have 
assumed that the precession frequency of the electron is precisely defined and there are no resonance offsets. The 
dispersion in the frequency of the electron can be refocused by the following techniques. For k even, let ri — vr — 
T2 — • • ■ — Tfc — TT, be a sequence that implements the propagator Uf — exp(— Then the sequence P — P, where 
P = Ti — TT — T2 — • • • — Tfe, will implement the propogator exp{—iblq) and refocus all resonance effects. 

We now consider the problem of synthesizing a general rotation of the kind exp{—i{aS zip)) , in the presence of 
resonance offsets. Let Iq be orthogonal to Ip. Let us consider a switching sequence that implements exp(— iTr/q) 
compensating for offset effects and a switching sequence (assuming no resonance offset), Q = ti — tt — r2 — • • • — t^, 
such that k is even and Q — n implements exp(— i^^z/p). Then the sequence Q ~ R — Q — R will implement the 
desired propagator exp(— iaS'^/p) and compensate for offset effects on the electron. 

Remark 4 (Polarization Transfer from Electron to Nucleus:) The above described methods can now be used 
to transfer polarization from electrons to the nucleus. The desired transformation 

can be performed as follows. 

cxp(-tf 5„)^ cxp[-i^2S^Iy) cxp(— if5^)^ cxp(-if 25^/x)^ , . 

where the various unitary transformations can be performed in a way that are robust to resonance offsets. Furthermore 
all Bilinear Hamiltonians above can be synthesized by simply doing tt pulses on the electron as described before. 

Remark 5 (Electron Coupled to Many Nuclei:) We now consider the coupling Hamiltonian of a single electron 
^, coupled to many spin i, nuclei. The Hamiltonian for the spin system takes the form [l| 

N N 
Hq = flsSz + ^ ^'jlljz + Sz ^(Aj/jz + Bjljx). 
3=1 J=l 

For now, we neglect the resonance offset fls, for the electron. As before, for each electron- nuclear coupling, we can 
define the unit directions {dja,djp) and the precession frequencies {ujja,i-jjf}). Then a straightforward computation 
shows that if for all j, the precession frequencies ujja 7^ i^jf} and the frequency pairs {ujja,LUjf}) ^ {ujka,'-^kfj) are 
distinct, then all the commutators generated by —iHi and ~iH2, where 

N N N N 

span the space {Ijp, Szljq}- Therefore any rotation on the individual nuclei is achievable by simply switching the 
electron state by tt pulses. These rotations combined with rotations on the electron are sufhcient to synthesize any 
unitary transformation on the coupled electron- nuclear spin system. Let Uaj and Ufsj represent the desired unitary 
transformation on the a and f3 manifold of the j'*'* nuclei. Let ri — tt — ■ • • — vr — Tfe — tt, (fc even) be the pulse sequence 
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that implements such a transformation. Then the switching times can be computed by iterative maximization over 
Ti ~ Tfci of a trace function hke 

^ ^r(C/^C/(/3,, T,) . . . U{a,,n)) + tr{U'p^U{a,,Tk) . . . U{(3,,n)), 
j 

where U{aj,T) represents rotation of the j^^ nuclei around the axis aj, with angular frequency ujaj for time r. 
Gradient algorithms for maximization of such trace functions were presented in (241] 



III. CONCLUSION 



In this article, we developed new techniques for control of coupled electron-nuclear spin system. We showed that 
the spin system at magnetic field strengths where the Larmor frequency of the nucleus is comparable to the hyperfine 
coupling strength can be controlled by simply performing a sequence of flips (each followed by a suitable delay) on the 
state of the electron spin. We showed that this mode of switched control, allows for universal control of the nuclear 
spin states. This switched control of nuclear spin states is much faster in settings where the Hyperfine coupling 
is much larger than the Rabi frequency of the nucleus. Combined with single spin rotations of the electron spin, 
the switched mode control can be used to synthesize any unitary transformation on the coupled spin system. We 
presented algorithms for explicitly computing the switching times. Furthermore, we discussed how these methods 
can be made robust to resonance offsets of the electron spin. Further method development is required to make 
these methods robust to dispersion in strength of hyperfine couplings. Application of these methods for transferring 
polarization from electron to nuclear spins was also discussed. The switched control methods presented in this work 
are expected to find applications in pulsed EPR experiments and quantum information processing systems based on 
coupled electron-nuclear spin dynamics. 
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